Self-organized quantum transitions in a spin-electron coupled system 
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We investigate quantum dynamics of the excited electronic states in the double-exchange model at 
half-filling by solving coupled equations for the quantum evolution of electrons and Landau-Lifshits- 
Gilbert equation for classical spins. The non-adiabatic quantum transitions driving the relaxation 
are coordinated through the self-organized space-time structure of the electron/spin dynamics lead- 
ing to a resonant precession analogous to the ESR process. 



In the strongly correlated electron systems such as 
the transition metal oxides, the rich phases including 
the spin/charge/orbital ordered states and supercon- 
ducting state are realized by the electron-electron and 
electron-lattice interactions An essential feature of 
the strongly correlated electron systems is the appear- 
ance of the internal degrees of freedom such as spins. 
The dynamics of the electrons with the strong interac- 
tion can be translated to the motion of the electrons 
in the background of the fluctuating spin fields via the 
Stratonivich-Hubbard transformation [l|. Therefore, the 
quantum dynamics of the correlated electrons can be for- 
mulated as that of the coupled dynamics of the spins and 
electrons. 

In particular, the time-dependent quantum dynamics 
in nonequilibrium states are now an issue of intensive 
interests. An example is the photo-induced quantum dy- 
namics of the spin-electron coupled systems 0, H, 0, H @ ■ 
In the systems, the dynamics of the spins plays a crucial 
role in the relaxation of the electrons after the photo- 
excitations leading to distinct features compared with 
the usual semiconductors or band insulators, where the 
electron-phonon interaction plays the major role in the 
relaxation. For example, the photo-induced ferromag- 
netism can leads to the metallic state in manganites due 
to enhanced kinetic energy in the double exchange in- 
teraction Q in sharp contrast to the self-trapping of 
the small polaron in the band insulators. This means 
that the coupled dynamics of the spins and electrons 
must be treated as the many-body collective phenomena 
in the correlated electron systems. The nonequilibrium 
quantum phenomena still remain challenges for theoret- 
ical study 0, S Q- In general, it is not easy to cap- 
ture the nature of the many-body electronic state, and 
our study below is complementary to the previous works 
giving a semiclassical picture for quantum dynamics for 
larger size/higher dimensional systems. 

In this paper, we study the relaxation dynamics of the 
excited states in the double-exchange model where 
the classical spins are coupled to the conduction elec- 
trons. This model has the advantage that the fully quan- 
tum dynamics of the electrons combined with the classi- 
cal motion of spins can be simulated, providing a clear 
physical picture. This way of describing the correlation 
effect is essentially justified for the magnetically ordered 
state where the (staggered) magnetization behaves as 
semi-classical object. We found that the self-organized 



space-time structure of the spins and electrons is formed 
through the quantum transition as shown below. 

We start with the Hamiltonian on the square lattice, 

H^-t ^ c\sCjs + h.c. - JhY^ c^Cis'ass' ■ S„ (1) 

<ij>^s iss' 

where <ij> denotes a nearest-neighbor pair, s and s' are 
indices for electron spin, respectively, and ass' is given 
by Pauli matrices. The local spins, S'i's, are taken to 
be classical vectors with magnitude S. Other notations 
are standard. We consider the half-filled case, i.e., one 
electron per site. 

Using finite size systems, we numerically investigate 
the time evolution of the electronic states and local 
spins. The equation of motion for the local spins is 
expressed by the Landau-Lifschitz-Gilbert (LLG) equa- 
tion. Si — —Jni'^i) X S'i + aSi x Si, where {ui) 
is the expectation value of electron spin at site i, 
and the Gilbert damping constant Sa describes all 
the other relaxation processes of the spins than the 
coupling to the conduction electrons. When {cfi) 
is fixed, the solution of the LLG equation is given 
by 0,(r) = -[{SJH/tW,)\/[l + {Saf)]{t/S)T, and 
ei{T) = 2tan-i[tan(6'io/2)exp(S'a(/),)], where {6,,(j).i) is 
the polar and azimuthal angles of Si in the local spin co- 
ordinate where z axis is parallel to {ui) and initial value 
of (pi is zero, and T is time. The initial value of 6i is 
written by OiQ. Evolution of the electronic state \^{T)) 
is given by |$(T)) = f7(T)|$(0)), where U{T) is a uni- 
tary operator for the time evolution. If {Si} are fixed, we 
have U{T) — exp{—iHT). We successively calculate the 
evolution of the electronic state and the LLG equation 
for a small time increment AT to investigate evolutions 
of the system in total. The changes in {Si} and \^{T)) 
are reflected to the calculation through H and 

The Hamiltonian ([T]) is expressed by a bilinear form 
of fermion operators. In such a case, it is known that 
the electronic state remains to be a single Slater 

determinant state if the initial state is so [ll| . This allows 
us to increase the system size, which is advantageous for 
investigation of low-frequency dynamics with accuracies. 

Figure [1] shows the calculated result on the system of 
size 8x8 with the periodic boundary condition. As a 
typical example, a parameter set, t—1, SJh='2, Sa=l, 
S=l, Ar=0.008, is used. The lower panel of Fig. [^a) 
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is the time (T) dependence of the energy level structure. 
The Fermi level is taken to be zero. At around r~0, we 
clearly see the energy gap 25' Jh between upper and lower 
energy bands in the lower panel of Fig.[lja). We prepare 
the initial state in the following way: The ground state 
of the double-exchange model ([T]) at half filling is the an- 
tiferromagnetic (AF) insulating state because of the per- 
fect nesting-condition in this system. In order to mimic 
the thermal fluctuation, we introduce a random tilting 
of each spin from the AF configuration up to 0.1 radian 
which corresponds to the state with an excitation energy 
of '^O.OOlt from the ground state. The energy band is di- 
vided into upper and lower ones. For the initial state, we 
fill up the lower energy band, and then transfer two elec- 
trons from the lowest eigenstates in the lower energy band 
to the two highest levels in the upper energy band. In the 
upper panel of Fig. [Ija), the number of electron of the 
highest (lowest) energy states in the upper energy band 
is shown by the broken (dotted) line, and the solid line is 
the number of electron in the upper energy band. Figures 
[IJb)-(e) show the time dependence of the configuration 
of the local spins, while Figs.[T](f)-(i) are the correspond- 
ing local energy density defined by the expectation value 

of -(i/2)Ep,s(4sCps + h.c.) - JnY^ss' ^sCis'^ss' ■ S^, 
where p runs over the nearest-neighbor sites of i. In the 
Figs. [IJf)-(i), the energy of the ground state is taken to 
be zero. 

In the numerical simulation, we find several distinct 
time regions. Stage (I)-(IV): 

Stage (I): Precursor process. Up to T'^20, there occurs 
almost no change and the spins are almost frozen at the 
initial configuration. Figure [TJf ) shows the energy den- 
sity at T=8, showing that the excitation energy is almost 
uniformly distributed due to the extended nature of the 
excited states. However, small deviation from the AF 
configuration grows eventually leading to the next stage. 

Stage (II): Self-organization process. At this stage, the 
system shows spatial inhomogeneity with energy dissi- 
pation. At around T=20~40, the spins start to move, 
and the electrons at the highest 2 energy states start to 
spread into lower energy states within the upper energy 
band (see the broken line in the upper panel of Fig.[lja)). 
Up to T^170, the deviation of the spins are not so large 
moving around the original direction as shown by the 
snapshot in Fig. [TJc) at T=80. This small amplitude 
spin fluctuation can induce the intra-band transitions of 
the excited electrons to lower and lower energy state. 
It is noted that the total energy decreases due to an 
energy dissipation originate from the Gilbert damping. 
The energy distribution shows the gentle spatial depen- 
dence with the reduced average as shown in Fig. [Ijg) for 
T = 80. This means that the electronic wavef unctions are 
rather extended though slightly disturbed by the small 
tilting of the spins. 

Stage (III): Relaxation process with interband transi- 
tion. The inhomogeneity developed in the previous stage 
brings about a remarkable localization behavior and de- 
rives a dynamic relaxation process with the interband 
transition. At around r~220, the electronic and local- 
spin structures show a dramatic change characterized by 
the large amplitude motion of a local spin marked by an 
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FIG. 1: (Color online) Time evolution of the double exchange 
model (see text), (a) The lower panel shows the time (T) 
dependence of the energy level structure. In the upper panel, 
the number of electrons in the upper energy band (solid line), 
in the highest two energy states (broken line), and in the 
lowest energy states of the upper energy band (dotted line) 
are shown. Configuration of local spins at r=8, 80, 224, and 
288 are presented in (b), (c), (d) and (e), respectively. Dot 
indicates the head of local spin. The local spins marked by 
ellipses in (d) and (e) are strongly deviated from the ground 
state. Figures (f), (g), (h) and (i) show the distribution of 
excitation energy density measured from the ground state at 
r=8, 80, 224, and 288, respectively. 



ellipse in Fig. [TJd). This motion starts at r-200. At 
this time, the excited electrons reach the lowest states 
in the upper energy band as shown by the dotted line in 
the upper panel of Fig. [T{a), and find the localized place 
to relax furthermore. That is, the localization of the 
electronic state together with the large-amplitude local- 
spin motion with the polar angle of the order of tt oc- 
curs concomitantly. This corresponds to a pair of the 
in-gap energy levels split off from the edge of the up- 
per and lower energy bands at around r=200, as seen 
in the lower panel of Fig. [ija). When the separation of 
in-gap energy levels is smallest, the number of excited 
electrons decreases rapidly. After that, the local spin re- 
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covers toward its original direction, and the in-gap states 
merges again into the upper and lower energy bands. As 
shown in Fig. [ijh) for T—224, the excitation energy is 
concentrated around this local spin, which in turn drives 
a motion of that local spin. In the same way, the relax- 
ation dynamics with interband transition occurs again at 
T~300, around another site (see Figs.[lja), (e) and (i)). 
This transition through the in-gap states reminds us the 
Landau-Zener mechanism. However, a more thorough 
study given below reveals that it is a resonant transi- 
tion and is completely different from the Landau-Zener 
process. 

Stage (IV): Relaxation process to a meta-stable state. 
After r^400, the active motion of the local spins is fin- 
ished and the alignment becomes nearly perfect AF. The 
excited electrons, however, remains more than ~0.7 in 
the upper energy band. This meta-stable state contin- 
ues for a long time within our simulation (at least up to 
T-8000). 

Now we consider the quantum dynamics in more depth. 
As discussed below, there are two components of the lo- 
cal spins, i.e., the rapid oscillation and the slow motion 
as expressed by Si—Sf°'^+Sl"'^^'^. Let us first consider 
the rapid oscillation. Figure[2^a) shows the y-component 
of the local spin moving with a largest deviation from 
the AF ground state configuration through Stage (I) and 
the early period of Stage (II), and the inset is the tra- 
jectory of the local spin on the Sx—Sy plane. As seen in 
Fig. [2l^a) , the local spin shows an oscillation with a period 
of Tp^l2, i.e., the frequency u}{=2it /Tp)^{).b. We find 
that this frequency w corresponds to the difference of the 
energy between the highest (ei) and second highest (£2) 
energy levels in the lower panel of Fig. [TJa). The rapid 
oscillation §1°'^^'^ is driven by time dependence of a in 
the LLG equation: When the wavefunction has the form 
\^{t))= ci(t)|l)+C2(i)|2) withc,(i) = Ca(0)e-*^"*, {a^it)) 
has the component proportional to ci{t)*C2{t){l\(Ji\2) 
cx e'*^^^"^^-** and its complex conjugate. Putting this 
into the LLG equation, we obtain Sl°'^^'^{t) cx Ci{t)*C2{t) 
(l|i7i|2) X -|- h.c This interpretation is consistent 

with the Fourier spectral weight '1^ of the spin motion 
in Fig. m^e) , where the peak is observed around the fre- 
quency Lo = 0.5, which corresponds to Si — £2 in Stage 
(I). This rapid oscillation in turn induces the transition 
between the state |1} and |2}, analogously to the elec- 
tron spin resonance (ESR) where the oscillating trans- 
verse magnetic field induces the Rabi oscillation. Figure 
[2jb) is the time dependence of the electron occupation 
number at the highest energy states at Stage (I)- (II), and 
it clearly shows this Rabi oscillation with the frequency 
(fi) determined by the oscillation amplitude (55""°^'*'^) of 
grapid ^ From Fig. Efa), we can read the oscillation am- 
plitude (55'''^^*''=0.3, so that the frequency is estimated 
to be JhSS^'^p^'^ =0.6. Therefore, the occupation num- 
ber will show an oscillation with a period of 27r/(217)=5.2. 
This oscillation is actually observed as shown in Fig.[2jb). 
With the Gilbert damping, this oscillation is the damped 
one and the occupation number of the lower energy state 
increases. 

It is important to note that the spatial inhomogeneity 
essentially occurs with the relaxation dynamics discussed 
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FIG. 2: Relaxation driven by rapid oscillation, (a) The y- 
component of the local spin which shows the largest deviation 
from the AF ground state configuration in Stage (1) and (11). 
The inset is the trajectory of the local spin on the Sx — Sy 
plane, (b) The electron occupation number at the highest 
energy states, (c) The same with (a) but in Stage (III), (d) 
The electron occupation number at the lowest energy state 
(in-gap state) in the upper energy band, (e) The Fourier 
spectral weight of (a), (f) The same with (e) but T—132 in 
Stage (11). (g) The Fourier spectral weight of (c). (h) The 
same with (f) but T=600 in Stage (IV). 



above. At Stage (I), the local spin alignment is almost 
perfect AF. In this state, the momentum is a good quan- 
tum number and the eigenvalue and eigenstate are well 
defined by it. With the time evolution, the excited elec- 
trons are distributed into lower energy levels, so that the 
spatial inhomogeneity appears as observed in the spatial 
distribution of the local energy density [13] ■ Refiecting 
the inhomogeneity, the magnitude of the rapid oscilla- 
tion of the local spin depends on the sites in real space 
strongly. When we look at the local spins at other sites, 
this rapid oscillation is almost missing and only the slow 
and small amplitude motion is observed. This means that 
the site-selective lock-in of the rapid oscillation of the 
spins occurs self-consistently with the electronic levels be- 
fore and after the transitions. This self-organized space- 
time structure is the most basic mechanism of the quan- 
tum transitions in the spin-electron coupled systems. 

With further time evolution in Stage (II), the en- 
ergy level structure becomes disordered reflecting the 
disordered spin configuration, and the spectral denstiy 
does not show characteristic frequency in this case (see 
Fig. m^f)) since it is given by many components corre- 
sponding to various e„ — £m- 

Figure [^I^c) shows the motion of the local spin marked 
by an ellipse in Fig. [TJd) for Stage (III). As seen in the 
inset of Fig.^^c), the local spin shows a rapid precession. 
Figure 2(d) shows the electron occupation number at the 
lowest energy state (in-gap state) in the upper energy 
band. In this case, however, the Rabi-oscillation behav- 
ior has not been observed. This is because the excited 
electrons in the states forming the bottom of the upper 
energy band show a cascade relaxation process, and the 
several frequencies are involved in those (see Fig.[5I^g)). It 
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FIG. 3: (Color online) Time evolution of the double exchange 
model. The parameter set, t=l, SJh='2, 5'q:=0.01, 5=1, 
Ar=0.008, is used, (a) The same with Fig. [TJa). (b) Struc- 
ture of local spins at r=2000. Distribution of energy density 
at (d) r=2000, (e) T=2100, and (f) T=2200. 

is seen, in fact, that at the early period of Stage (HI), the 
excited electrons occupy the lowest energy states of the 
upper energy band, and from the those states an in-gap 
state appears (see Fig. [1]) . Corresponding to the differ- 
ences between those energy levels, the spectral density in 
Fig. dig) has the peak around u = 1. In other words, the 
inter-band transition occurs successively through those 
energy levels. Therefore, the dynamics similar to the 
ESR process is also essential for the inter-band transi- 
tion in this Stage (III). In Stage (IV), there occurs no 
quantum transition any more because (l|(Ti|2) xSf°'^=0 
for the AF state. 

Let us consider the time scale for the slow motion of 
the local spins. The deviation of the electron spin from 
the ground state of the excited electronic state gives a 
force to the local spins through the LLG equation. There- 
fore, {(l/7Veff) Jfl-S'a/[1 + (Sa)^]}^^ determines the time 
scale for the motion of the local spins, where A'cfi is the 
number of the sites in the real space over which the wave- 



function of the excitation has the finite weight (see the 
solution of the LLG equation) . For the plane wave states 
at Stage (I) and (II) , Nctf = N (the total number of the 
sites) and decreases as the wavefunction becomes more 
localized with the time evolution as seen in Figs.[T][f)-(i). 
In the present case, Sa = 1, a time scale of the slow dy- 
namics being of the order of 100 is derived We have 
also studied the more realistic case of S'a=0.01 (Fig. [3]). 
The mechanism of the relaxation dynamics for 5*0=0.01 
is essentially the same with the previous case. The ob- 
served time scale of the slow dynamics in Fig. [3]Ja) is 
about 10 times as large as that of Fig. [Ija), whereas the 
above simple consideration on the time scale gives about 
50 times longer one compared with the case that Sa — 1 . 
In the case for small Sa, because the friction is small, 
a large number of local spins show very active dynamics 
(see Fig.[3ljb)). This results that the considerable spatial 
area (but not the whole area of the system) and many 
in-gap states are responsible for the relaxation dynamics 
in Stage (III) (see Figs. [31Ja), (d) and (e)), whereas the 
strongly localized in the real space and a few in-gap states 
are available for the relaxation dynamics in the case that 
Sa = 1 (see Figs. H^a), (h) and (i)). As a result, the 
relaxation process is accelerated, and the 10 times differ- 
ence in the slow time scale appears although a hundred 
times difference in the magnitude of Sa is in this case. 

In summary, we have studied the relaxation dynamics 
of the excited state of the electron-spin coupled systems. 
The self-organized space-time structure of the spins and 
electronic state triggers the quantum transitions, and the 
adiabatic approximation does not work here. 

The authors are grateful to Y. Tokura, M. Kawasaki, 
H. Matsueda, T. Tohyama, S. Ishihara, and K. Tsutsui 
for useful discussions. This work is supported by Grant- 
in-Aid for Scientific Research (Grant No. 19048015, 
19048008, 17105002, 21244053, 18560043, 19014017, and 
21360043) and a High-Tech Research Center project for 
private universities from the Ministry of Education, Cul- 
ture, Sports, Science and Technology of Japan, Next 
Generation Supercomputing Project of Nanoscience Pro- 
gram, JST-CREST and NEDO. 



[1] M.Imada, A. Fujimori, and Y. Tokura,, Rev. Mod. 

Phys.70, 1039 (1998). 
[2] Y. Ogawa, S. Koshihara, K. Koshino, T. Ogawa, C. 

Urano, and H. Takagi, Phys. Rev. Lett. 84, 3181 (2000). 
[3] Y. Mitsumori, A. Oiwa, T. Slupinski, H. Maruki, Y. 

Kashimura, F. Minami, and H. Munekata, Phys. Rev. 

B 69, 033203 (2004). 
[4] S. Koshihara, A. Oiwa, M. Hirasawa, S. Katsumoto, Y. 

lye, C. Urano, and H. Takagi, and H. Munekata, Phys. 

Rev. Lett. 78, 4617 (1997). 
[5] K. Miyano, T. Tanaka, Y. Tomioka, and Y. Tokura, Phys. 

Rev. Lett. 78, 4257 (1997). 
[6] M. Fiebig, K. Miyano, Y. Tomioka, and Y. Tokura, Sci- 
ence 280, 1925 (1998). 
[7] H. Matsueda and S. Ishihara, J. Phys. Soc. Jpn. 76, 

083703 (2007). 

[8] H. Matsueda, A. Ando, T. Tohyama, and S. Maekawa, 
Phys. Rev. B 77, 193112 (2008)'. 



[9] K. Yonemitsu and N. Maeshima, Phys. Rev. B 79, 125118 
(2009). 

[10] P.-G. de Gennes, Phys. Rev. 118, 141 (1960). 
[11] M. Imada and Y. Hatsugai, J. Phys. Soc. Jpn. 58, 3752 
(1989). 

[12] For the Fourier analysis, 2^^ points with a small time 
increment AT^O.OOS in units of 1/t are used. 

[13] We have examined how long it takes to begin the inter- 
band transition in Stage (III) using 6x6, 8x8, 10x10 and 
12x12 systems. In the initial state, two electrons are ex- 
cited from the lowest to highest energy levels. We have 
confirmed that the period leading up to the inter-band 
transition increases linearly as a function of the system 
size. 

[14] The inhomogeneity has also been observed in the spatial 
distribution of the magnitude of electron spins 



